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Superfluid states of symmetric nuclear matter with finite total momentum of Cooper pairs (nuclear 
LOFF phase) are studied with the use of Fermi-liquid theory in the model with Skyrme effective 
forces. It is considered the case of four-fold splitting of the excitation spectrum due to finite 
superfluid momentum and coupling of T = and T = 1 pairing channels. It has been shown that at 
zero temperature the energy gap in triplet-singlet (TS) pairing channel (in spin and isospin spaces) 
for the SkM* force demonstrates double-valued behavior as a function of superfluid momentum. As 
a consequence, the phase transition at the critical superfluid momentum from the LOFF phase to the 
normal state will be of a first order. Behavior of the energy gap as a function of density for TS pairing 
channel under increase of superfluid momentum changes from one-valued to universal two-valued. 
It is shown that two-gap solutions, describing superposition of states with singlet-triplet (ST) and 
TS pairing of nucleons appear as a result of branching from one-gap ST solution. Comparison of the 
free energies shows that the state with TS pairing of nucleons is thermodynamically most preferable. 
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In this report we shall study superfluid states of nu- 
clear matter with nonzero total momentum of Cooper 
pairs. First the states with moving condensate were con- 
sidered in Rcfs. jlj H with respect to metallic supercon- 
ductors. In this case the superconducting condensate 
has a spatially periodic structure. Corresponding phase 
is called the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) 
phase. Recent upsurge of interest to the LOFF phase is 
related with the possibility of formation of this phase in 
nuclear Q and quark matter (J, || . In Rcf . || it was con- 
sidered in model calculations the case of neutron-proton 
superfluidity, when the quasiparticle spectrum is two- 
fold split due to asymmetry of nuclear matter and finite 
superfluid momentum. It was shown, that the nuclear 
LOFF phase appears as a result of a first order phase 
transition in the asymmetry parameter from the spatially 
uniform BCS superconducting state. Further increase of 
the asymmetry parameter leads to the second order phase 
transition from the LOFF phase to the normal state. In 
the study of Ref. || it was assumed that coupling be- 
tween isospin singlet and isospin triplet pairing channels 
can be neglected. However, as emphasized in Ref. 0, in 
the region of low densities coupling between T = and 
T = 1 pairing channels may be of importance, leading 
to the emergence of multi-gap superfluid states, charac- 
terized by nonvanishing gaps in both pairing channels. 
Thus, we shall consider the case of four-fold splitting of 
the quasiparticle spectrum, caused by account of, first, fi- 
nite superfluid momentum and, second, coupling of T = 
and T = 1 pairing channels. Another simplifying mo- 
ment in Ref. || is the use of free single particle spectrum 
and bare interaction in the gap equation without taking 
into account the effects of medium polarization. A "first- 
principle" derivation of the pairing interaction from the 
bare NN force still encounters many problems, such as, 
e.g., treatment of core polarization W. Hence, it is quite 
natural to develop some kind of a phenomenological the- 



ory, where a phenomenological pairing interaction is em- 
ployed. As such a theory, we shall use the Fermi-liquid 
(FL) approach ||. In the Fermi-liquid model the normal 
and anomalous FL interaction amplitudes are taken into 
account on an equal footing || . This will allow us to con- 
sider consistently influence of the FL amplitudes on su- 
perfluid properties of the nuclear LOFF phase. Besides, 
as a potential of NN interaction we choose the density de- 
pendent Skyrme effective forces, used earlier in a number 
of contexts for description of superfluid properties both 
finite nuclei jlO], and infinite nuclear matter |l2|, . 

Superfluid states of nuclear matter are described 
by the normal f KlK2 = Tr ga+ 2 a Kl and anomalous 



Tr ga K 



distribution functions of nucleons 



(k = (p, er, t), p is momentum, <j(t) is the projection 
of spin (isospin) on the third axis, g is the density matrix 
of the system) . We shall study two-gap superfluid states 
in symmetric nuclear matter, corresponding to superpo- 
sition of states with total spin S and isospin T of a pair 
S = 1, T = (triplet-singlet (TS) pairing) and 5 = 0, 
T = 1 (singlet-triplet (ST) pairing) with the projections 
S z = T z = (TS-ST states). In this case, assuming that 
a condensate moves with the finite momentum q, the nor- 
mal / and anomalous g distribution functions read |6| |^] 

/«i« 2 = foo(pi)(0oTo)KiK 2 Spi,P2> (1) 

.9kik 2 = (.93o(Pl)030-2T2 + ff03 (Pi )o2T 3 T 2 ) Kl K2 <5 pi ._p 2+g , 

(2) 

where <Ji and are the Pauli matrices in spin and isospin 
spaces, respectively. For the energy functional, being 
invariant with respect to rotations in spin and isospin 
spaces, the quasiparticle energy and the order parameter 
have the similar structure 
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If to take into account the antisymmetry properties A T = 
—A, g T = —g and to set Pi = p + q/2, P2 = p + q/2 
(q is total momentum of a pair, p is momentum of one 
of nuclcons in the center of mass frame of a pair), then 
we obtain 



A 30 (p + -j) = A 30 (-p 
A 03 (P + 7r) = A 03 (-p 



) = A 30 (p,q) 
) = A 03 (p,q) 



(5) 



and analogous relationships hold for the functions 
530j5o3- Further we shall write the self-consistent equa- 
tions for the quantities A 30 (p,q), A 03 (p,q). Using the 
minimum principle of the thermodynamic potential and 
procedure of block diagonalization [||], one can express 
evidently the distribution functions /ooj<?30i3o3 m terms 
of the quantities e and A: 
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with TS and ST pairings of nucleons the energy func- 
tional is characterized by one normal Uq and two anoma- 
lous Vi,V 2 FL amplitudes [| ||. Then one can obtain 
the self-consistent equations in the form 
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where mo being the mass of a bare nucleon. 

Further for obtaining numerical results we shall use the 
Skyrme effective interaction. In the case of Skyrme forces 
the normal and anomalous FL amplitudes read |p_3|| 



UqQl) = 6t + 6t 3 g 
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/i° is chemical potential, which should be determined 
from the normalization condition 



where U,Xi,(3 are some phenomenological constants, 
characterizing the given parametrization of Skyrme 
forces (we shall use the SkP jL0j, SkM* (l4| potentials). 
With account of Eq. (|T3|) we obtain 
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g is density of symmetric nuclear matter. As follows from 
the structure of the distribution functions /oo, <730> 303, 
the quantity 



and the effective chemical potential fj, should be deter- 
mined from the normalization condition (||). Besides, 
expression for the quantity £ a reads 



£«(±p+~) 



E±, 



being the exponent in Fermi distribution functions 
n±(±p + §), plays the role of the quasiparticle exci- 
tation spectrum. In the considering case the spectrum 
is four-fold split due to 1) finite superfluid momentum 
(q 0)) 2) coupling of TS and ST pairing channels 
(A 30 ^0,A 3^0). 

To obtain the closed system of equations for the quan- 
tities A and £, it is necessary to set the energy functional 
of the system. In the case of symmetric nuclear matter 
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The normal distribution function /oo in turn depends on 
the quantity £ a and, hence, expression ( [Tt| ) represents 
an equation for determining the quantity £ a . Since the 
second term in Eq. (17) is proportional to the scalar prod- 
uct pq, solution of Eq. (17) should be found in the form 



£a(p + f ) = Pq/2m*, where m* is some effective mass. 
Using Eqs. (^)-(^), we present equations for the order 
parameters A 3 q,Ao 3 , effective chemical potential [i and 
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effective mass m* as 
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Eqs. (|i"g|)-(|2l"|) describe two-gap superfluid states of 
symmetric nuclear matter with moving condensate and 
contain one-gap solutions with A30 7^ 0, A03 = (TS 
pairing) and A 30 = 0, A 3 7^ (ST pairing) as some 
particular cases. We shall analyze Eqs. (|l8|)-(^ll) using 
the simplifying assumption, that FL amplitudes V\ , V2 
are not equal to zero only in a narrow layer near the Fermi 
surface: \£ s \ < 6, 9 <C Sf (further we set 9 = O.lep). 

First we shall find the dependence of the order pa- 
rameters A3o(p = pf), Ao3(p = pf) from the superfluid 
momentum at zero temperature. We begin our analy- 
sis with finding one-gap solutions of the self-consistent 
equations. Results of numerical determination of the en- 
ergy gap are presented in the Fig. [l]a. It is seen, that 
a general tendency is quite clear: at low superfluid mo- 
menta the energy gap retains its constant value and then 
rapidly decreases and vanishes at some critical point. 
This means, that one-gap superfluid states with moving 
condensate will disappear at large enough superfluid mo- 
menta. However, for TS pairing of nucleons, interacting 
via SkM* effective potential, the phase curve has an inter- 
esting peculiarity, namely, the energy gap demonstrates 
nontrivial double-valued behavior in the region, where 
it sharply falls. Such behavior differs from the ordinary 
one-valued behavior of the energy gap, as, e.g., in the 
case of TS pairing and SkP effective interaction. To un- 
derstand this difference, we have determined the effective 
mass m* as a function of superfluid momentum, Fig. |]b. 
For the SkP potential the mass m* is equal to the bare 
mass mo of a nucleon practically for all superfluid mo- 
menta. Unlike to this behavior, the effective mass m* 
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FIG. 1: Energy gap (top) and effective mass m* (bottom) 
vs. superfluid momentum for different types of pairing and 
Skyrme forces at g = 0.03 fm -3 and zero temperature. 



for the SkM* potential rapidly decreases close to the re- 
gion of phase transition to the normal state. Since the 
mass m* enters into Eqs. (p^|)-(^ll) only through the ratio 
q/m*, descent of the effective mass leads to the increase 
of the effective shift between the centers of proton and 
neutron Fermi spheres, on which the paired proton and 
neutron lie. This gives the possibility to the appearance 
of the second solution for the energy gap. According to 
Eq. (|l5|), the parameter eo of the normal FL amplitude 
Uq determines the sign before the sum in the l.h.s. of 
Eq. (pl|), and, hence, determines whether the mass m* 
will be greater (if eo < 0) or less (if eo > 0), than the 
bare mass mo (if corresponding sum is nonzero). For 
the SkM* potential e > while for the SkP potential 
eo < 0, that explains the difference in the behavior of the 
effective mass for these two forces. However, diminishing 
behavior of the effective mass m* does not guarantee that 
the energy gap will have two-valued behavior as a func- 
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tion of superfiuid momentum. In Fig. la we have plotted 
also the dependence of the energy gap for ST pairing of 
nucleons and SkM* potential. As follows from here, in 
spite of correctness of the inequality eo > 0, the energy 
gap demonstrates the usual one-valued behavior. Since 
ST coupling constant in the Skyrme model is always less 
then TS one ||, one can conclude, that the pairing inter- 
action should be strong enough for the second solution 
to be developed. The possible two-valued behavior of 
the energy gap has important consequence. If to go from 
the region of large enough superfiuid momenta in the di- 
rection of low momenta, then at some critical point the 
energy gap will arise by a jump, and, hence, the phase 
transition to the superfiuid phase will be of a first order. 

Now we go to the study of two-gap superfiuid states 
when both order parameters, A 30 and A 03 , are not equal 
to zero. Results of numerical determination of the or- 
der parameters A 30 (q), Ao 3 (q) on the base of Eqs. (|i~8|)- 
( pl| ) are presented in Fig. |[ Here tsst(ts) and tsst(st) 
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FIG. 2: Order parameters A30, A03 vs. superfiuid momentum 
for SkM* force at g = 0.03 fm -3 and zero temperature. 

are notations for the dependencies of TS and ST order 
parameters in the TS-ST solution of the self-consistent 
equations. As seen from Fig. ||, TS-ST solutions appear 
as a result of branching from one gap ST solution (in 
the branching point A30 = 0, A03 = Aq 3 , Aq 3 being one- 
gap ST solution). Note that the self-consistent equations 
have two-gap solutions in the case of SkM* potential, but 
have no such solutions for the SkP potential. As clarified 
in the Ref . [|| , for the existence of TS-ST solutions it is 
necessary that TS and ST coupling constants must be of 
the same order of magnitude. However, this condition is 
broken for the SkP potential, where TS coupling constant 
is much larger than ST one. 

Since we use density dependent effective interaction, it 
allows us to study also the dependence of the order pa- 
rameters from density of nuclear matter. The results of 
numerical determination of one-gap and two-gap solu- 



tions of self-consistent equations at fixed superfiuid mo- 
mentum are presented in Fig. 0. As one can see, super- 



SkM" 




0.14 



FIG. 3: Order parameters A30, A03 vs. density for SkM* force 
at q/pF — 0.04 and zero temperature. 

fluidity with finite superfiuid momentum exists in finite 
density region (gi(q), ^2(9)), excluding some vicinity of 
the point g = (for TS pairing the left point g± (q) is very 
close to g = 0). The most important peculiarity, i.e., the 
double-valued behavior of the energy gap for TS pairing 
in the case of SkM* potential is preserved for the given 
ratio q/pF- For other types of pairing and Skyrme forces, 
considered earlier (including the SkP potential, which is 
not shown in Fig. |3|), the order parameter for one-gap so- 
lutions has one-valued behavior. In the case of two-gap 
solutions the mechanism of their appearance is similar to 
the considered above, i.e., it is branching from one-gap 
ST solution. 

However, in general case behavior of the energy gap 
as a function of density is more complicated. In Fig. |4J 
we plot the dependence of the energy gap in the case of 
TS pairing and SkM* interaction for the set of fixed val- 
ues of the ratio q/pF- It is seen, that behavior of the 
phase curves may be one-valued or two-valued, that de- 
pends on the value of the ratio q/pF- At small enough 
superfiuid momentum (e.g., for q = 0.009p/) the gap be- 
haves as a one-valued function of density; for the ratios 
q/pF, larger than some critical value qi/pF, the gap has 
two-valued behavior in the region close to the right crit- 
ical point g2{q)- Further increase of the ratio q/pF leads 
to formation of the second part of the phase curve with 
double-valued behavior in the region close to the left crit- 
ical point gi(q) (e.g., for q = O.lpf). When q increases, 
the regions with double-valued behavior of the gap begin 
to approach and at q — q c (q c « 0.106j?f) it takes place 
contiguity of the regions with two solutions. For q > q c 
the phase curves are separated from the density axis and 
turn into the closed oval curves. Under further increase of 
q the dimension of the oval curves is reduced and at some 
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FIG. 4: Energy gap vs. density in the case of TS pairing and 
SkM* force for different values of the ratio q/pF- 

1 = 1m the oval curves shrink to a point (q m w 0.108pf)- 
For the values q > q m TS superfluidity of nuclear mat- 
ter vanishes. Thus, in the range < q < q\ the gap is 
one-valued function of density, in the range q\ < q < qi 
the phase curve A30 = A3o(g) has one part with double- 
valued behavior, for q-i < q < q c it contains two distinct 
parts with double-valued behavior and for q c < q < q m 
the gap has a universal double- valued behavior. 

Since we have a few solutions of self-consistent equa- 
tions, it is necessary to check, which solution is thermo- 
dynamically favorable. Calculations show that due to 
the large size of the gap in TS pairing channel the free 
energy of the corresponding state as a function of momen- 
tum or density much smaller then for the case of ST and 
TS-ST pairing. Hence, TS superfluid state wins compe- 
tition for the thermodynamic stability. If to compare the 
free energies of two different branches, corresponding to 
double-valued behavior of the energy gap in TS pairing 
channel, then the branch with larger size of the gap will 



be thermodynamically favorable. 

In summary, we studied superfluidity of symmet- 
ric nuclear matter with moving condensate in the FL 
model with density dependent Skyrme effective interac- 
tion (SkM*, SkP forces). It has been considered the case, 
when the quasiparticle excitation spectrum is four-fold 
split due to finite superfluid momentum and coupling of 
T = and T = 1 pairing channels. Apart from the 
renormalization of the chemical potential and bare mass 
of a nucleon, taking into account of the normal FL am- 
plitude leads to appearance of additional effective mass 
m* in the linear on superfluid momentum term in sin- 
gle particle energy. It is shown that at zero temperature 
the energy gap in TS pairing channel for SkM* poten- 
tial demonstrates two-valued behavior as a function of 
superfluid momentum. This is caused by the decreasing 
behavior of the effective mass m* close to the region of the 
phase transition to the normal state and strong enough 
interaction in TS pairing channel. The behavior of the 
energy gap as a function of density in TS pairing channel 
in general case is more complicated and under increase 
of superfluid momentum it changes from one-valued to 
universal two-valued character (until superfluidity disap- 
pears at some critical momentum). Two-gap solutions of 
self-consistent equations, corresponding to the case when 
both TS and ST order parameters are not equal to zero, 
appear as a result of branching from one-gap ST solution. 

Calculation of the free energy shows that TS superfluid 
state is thermodynamically most preferable state. In the 
case of double-valued behavior the gap changes in the 
critical point by a jump and, hence, the phase transition 
from the LOFF phase to the normal state will be of a first 
order (in superfluid momentum or density). Since the 
possible two-valued behavior of the gap will be preserved 
for small asymmetry, this will be also true for weakly 
asymmetric nuclear matter, that differs from the picture 
of a second order phase transition, considered in H. 
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